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Abstract 

We establish a fundamental connection between smooth and polygonal knot ener- 
gies, showing that the Minimum Distance Energy for polygons inscribed in a smooth 
knot converges to the Mobius Energy of the smooth knot as the polygons converge 
to the smooth knot. However, the polygons must converge in a "nice" way, and 
the energies must be correctly regularized. We determine an explicit error bound 
between the energies in terms of the number of the edges of the polygon and the 
Ropelength of the smooth curve. 
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1 Introduction 



Given a knot K in 3-space, there are several ways to define an "energy func- 
tion" that measures how complicated the knot is in its spatial conformation. 
In this paper, we establish a fundamental approximation theorem, showing 
that when both are appropriately normalized, the Minimum Distance Energy 
for polygons inscribed in a smooth curve converge to the Mobius Energy of the 
curve as the polygons converge to the smooth knot. We do a careful analysis, 
and determine an explicit error bound (Theorem 1), from which the approxi- 
mation (Theorem 2) follows immediately. 

In Section 2, we state the main theorems and agree on notation for the whole 
paper. In Section 3 we present a number of lemmas. These establish useful 
properties of curves and chords, so they may be of independent interest. In 
Section 4, we outline the proof of the error bound (Theorem 1), especially how 
to divide the problem into several cases (more precisely, divide the domains 
into different "zones") for which different analyses are needed. In Section 5, we 
give the detailed analyses for the various cases, and in Section 6, we combine 
the results from Section 5 to obtain the overall bound. 

Of course the error depends on how well the polygon approximates the smooth 
curve. However, there are more subtle issues to confront in controlling the 
error: One must reckon with the amount of curvature the knot has, and how 
close it is to being self-intersecting. These are captured by the thickness radius 
r(K) (see later in this section for definition). Our error bound is developed 
in terms of the total arc-length £(K), the number of edges of the inscribed 
polygon n, the mesh size 5 = the thickness radius r(K), and the ratio 
E L {K) = Since these quantities are interrelated, there are various ways 
to write the bound: the one we give in Theorem 1 is stated in terms of n and 
El(K) to emphasize that it is invariant under change of scale. 

Let t — > x(t) be a unit-speed parameterization of K with domain a circle C. 
The Mobius Energy or O'Hara Energy is 

E ° iK) = JLc i» w -irhp dsdL 

The energy E was defined and studied in jl|; 0: 0; 0) • The subscript in E 
reminds us that this version of the energy is exactly zero if K is a circle. 

By visualizing a smooth knot as being made of some "rope" , with a positive 
thickness, we obtain a fundamental measure of knot complexity. Hold the 
core knot K fixed and thicken the rope until the moment of self-contact. 
Call that sup radius the thickness radius or injectivity radius, r(K). Here is 
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a more precise definition: For small enough r, the knot K has a solid torus 
neighborhood consisting of pairwise disjoint disks of radius r centered at the 
points of K and orthogonal to K at those centers. Gradually increase r until 
some meridional disks touch; we call that supremum of good radii r(K). The 
ratio 

. arc-length of K 

El{K) = HK) ' 

called the Rope-Length of K, is a scale-invariant numerical measure of knot 
compaction. 



The basic theorems on thickness appear in (J5|), although the energy El first 
appeared in (Jy). We recall the properties of El in Section 3.3. 

Let P be a polygon with n edges. The Minimum Distance Energy of P is de- 
fined (0) as follows: For each pair X, Y of nonconsecutive edges of K, compute 
the minimum distance between the segments MD(X, Y), define U m d{X, Y) = 

lcngth(X) -length^) , 
[MD(X,Y)] 2 ' 

u' md (P)= E E u md (x,Y) . 

all edges X Yj^X or adjacent 

This version of U m d counts each edge-pair twice, analogous to a double integral 
over (most of) K x K. We write U' md just to distinguish from the original ver- 
sion (|7j) that counted each pair once. To consider knots with varying numbers 
of segments, we regularize by subtracting the energy associated to a standard 
regular n-gon (0) (or see (0)). Note that U' md is scale invariant, so we can use 
any regular n-gon and get the same number. We define 

E m d(P) = U' md (P) - C/^ (regular n-gon) . 



The energy U m< i has been implemented in several software systems (flCl : 11 ; 1^) , 
and studied in (Til 14; lih . 



We shall show that for suitable polygonal approximations P of a smooth curve 
K, E md (P) ps E (K). While the Mobius Energy is defined for C 1,1 curves, our 
proof requires that the knot K be C 2 smooth. 



In order for the approximation to work, we need to be careful about what it 
means to say, "the polygon P is a close approximation of K" . First, we need 
to prevent extreme changes in the edge lengths of P (see Figure 1). Suppose 
P is a polygon closely inscribed in K. We can slide vertex t> 3 along K towards 
vertex v^, making edge e% arbitrarily short, making edge C2 longer, and keeping 
the other edges of P fixed. This will make the contribution of the edge-pair 
(^2, e$) to E m d arbitrarily large. Thus, we can make polygons P' that also seem 
like close approximations of K, yet E md (P') » E (K). 
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Fig. 1. Extreme differences in edge lengths causes E m ^ to become large despite 
having the polygon close to the smooth curve. In such a case, E rn( i will be much 
larger than Eq. 



Fig. 2. Portions of the smooth curve can be arbitrarily close, causing Eq to be very 
large while E m( i of the inscribed polygon remains fixed. In such a case, Eq is much 
larger than E mc i- 

To prevent this problem, we have to limit the variation in edge lengths of 
polygons inscribed in K; we do this by having the vertices equally spaced in 
arc-length along K. One can modify our arguments to handle other tractable 
approximating polygons, e.g. equal edge lengths or "equal time" subdivisions 
of a regularly parameterized curve. 

Conversely, we can find situations where E m d(P) « E (K). Let P be the 
polygon (not drawn) (v±, v 2, t>3, V4, Vi) in Figure 2. We construct the quadrilat- 
eral so that the arcs between consecutive vertices are of equal length. Keeping 
the vertices fixed, deform the arc vyv-i and the arc V3V4 slightly so they get ar- 
bitrarily close to intersecting (where one crosses over the other in the figure). 
This makes E (K) — E rn( i{P) arbitrarily large. This problem is detected by 
the fact that r(K) decreases to 0, since normal disks of smaller and smaller 
radii will intersect. This is why the error bound in Theorem 1 must take into 
account the geometric quantity r(K). 

To avoid situations as in Figure 3, we assume the phrase "inscribed polygon" 
means the vertices of P occur in the same order as they occur along K. 

Finally, note that the regularizations play an essential role, making the proof 
more delicate than may be at first evident. See the discussion in Section 4. 
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Fig. 3. We require that the order of the vertices of the inscribed polygon coincides 
with an orientation of the smooth knot. 

2 Statement of main results and notation 



We shall use the following notation throughout the paper. 

K is a C 2 smooth simple closed curve in M 3 . 
£(K) is the total arc-length of K. 
r(K) is the thickness radius of K. 
El{K) = is the Rope-Length. 

5 = is the mesh size of the inscribed polygon. 
K is subdivided into n arcs of equal length. 5 = ^p-- 
Vi, . . . ,v n are the subdivision points along K. 

P n is the polygon formed by connecting the points v i, . . . , v n , v i in order. 
arc(x, y) is the length of the shorter of the two arcs of K connecting x and 

y- 

• \e\ is the length of the line segment e. 

Additional notation used in the proofs is listed at the beginning of Section 4. 

Theorem 1 (Error Bound) For any smooth knot K , if P n are inscribed 
polygons as above and n is large enough that n > E L (K), i.e. 5 < r(K), then 

\E (K)-E md (P n )\<$(n,E L (K)) , 

where $ is a linear combination (see final page of the paper) of six fractions 
of the form for various a > and b > 0. By combining some terms, 

we can take $ = 550 ^gf /4 + 10^r^. For very large n (n > E L (K) n / 3 ), 

we can use $ = 560 El< ^/\ 1 ■ 



Remark. There are other ways to write this scale-invariant error bound, using 
the identity 

E L (K) _ 5 
n r{K) ' 

From Theorem 1, we have immediately: 

Theorem 2 (Approximation Theorem) For any smooth knot K, if P n 
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are inscribed polygons as above, then as n — > oo ; E m d(P n ) — > Eq(K). 

PROOF. The supporting lemmas and the proof of the theorem occupy the 
rest of the paper. The lemmas are in Section 3. In Section 4, we outline the 
proof and explain how the domains will be divided into zones for which differ- 
ent analyses are needed. We give the analysis for each zone in Section 5 and 
put them all together in Section 6. 

There are numerous coefficients in the calculations; we constantly round up 
and pick the worst-case values, to keep the claims accurate and the numbers 
simple. 



3 The lemmas 

In this section, we prove the lemmas needed for the proof of the main theorem. 
3. 1 Lemmas about the cosine function, also chords and arcs of circles 
Lemma 3 If < <fi < n, then the following hold: 



PROOF. For (a), consider the Taylor series for cos(0). Parts (b), (c), (d) 
follow immediately. 

Lemma 4 

(a) On the unit circle C , for any points x,y, 



(a) 1 - \<f) 2 < cos0 < 1 - ±0 2 + , 

(b) 2 -^0 4 <2-2cos0<0 2 , 

(c) 1-^0 2 <^^<1, 

(d) ^<1 + |0 2 - 
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12 \x — y\ 2 &rc(x,y) 2 4 



(b) On a circle of radius R, 



111 1 




1 




12 R 2 \x - y\ 2 &rc(x,y) 



7T 



2 J R 2 ' 
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PROOF. Let be the angle (< ir) between points x and y on the circle. Since 
C is the unit circle, arc(:r, y) = and \x — y\ 2 = 2 — 2 cos0. The function 

1 1 

2 - 2COS0 ~~ 02 

is monotone, has a maximum at = tt, and is bounded below by the limiting 
value as — > 0. Part (b) is similar. 



Next we want to compare the quantities j^j^p an d jxZyp) where the points 
lie on circles of different sizes. 

Lemma 5 Suppose r < R are radii of circles and < a < nr. Construct any 
arcs of ( the same ) length a on the two circles and let x, y and X, Y be the 
endpoints of the two arcs. Then 

< \x-y\ 2 ~ \X-Y\ 2 < \A~^)r~ 2 ' 



PROOF. Chord length is always less than arc-length. For a fixed arc-length, 
as the radius gets larger, the chord length gets closer to the arc-length. Thus 
\X — Y\ > \x — y\. On the other hand, applying Lemma 4(b) to each circle, 
we have 

f 1 H 1 \ 1 1 1 

\x-y\ 2 ~ \X-Y\ 2 < \A~^ 2 ) ^'UB 2 ' 



3.2 Lemmas about chords and arcs of general curves 



We rely a lot on Schur's Theorem. Here is the version we need: 

Lemma 6 Let K be a C 2 smooth curve in M, 3 whose curvature everywhere 
is < some number k. Let C be a circle of curvature k, i.e. of radius r — K 
Let x,y G K,s,t G C such that aic(x,y) = arc(s,t) < nr. Then the chord 
distances satisfy 

\x — y\ > \s — t\ . 

When we write the chord length on C in terms of the central angle, this becomes 

( /arc(M)\\ 1/2 
\x-y\>r 2 - 2 cos 
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PROOF. See Schur's Theorem in (fl6h. 



Lemma 7 Let K be a C 2 smooth curve in M. 3 , with minimum radius of cur- 
vature r. Suppose x : [0, 7rr] — > M. 3 is a unit speed parameterization of an arc 
of K of length nr. Then the function \x(t) — x(0)\ is monotone increasing. 
That is to say: As points move farther apart along the curve, they also move 
farther apart in space, so long as the arc-distance is no greater than nr. 



PROOF. Let f(t) = \x(t) - x(0)\ 2 = (x(t) - x(0)) ■ (x(t) - x(0)). We claim 
§ > for t E (0, irr). The derivative f = 2(x(t) - x(0)) ■ x'(t). Thus we need 
to show that this dot product is positive, for all points x(t) in the interior of 
the arc. The proof uses the same central idea as the proof of Schur's theorem. 

We have 

x(t) — x(0) = / x'(s) ds , 
Jo 

so 

(x(t) — x(0)) ■ x'(t) = / x'(s) ■ x'Ct) ds . 

Jo 

The dot product x'(s) ■ x'{t) is just the cosine of the angle < n between the 
two velocity vectors. This angle is measured by the length of the geodesic 
arc on the unit sphere between the unit vectors x'{s) and x'{t). The trace of 
x'(u), as u runs from s to t, is another path on the unit sphere between the 
same vectors. The length of that path gives an upper bound for the length of 
the geodesic path. Thus, since < 1/r (recall r = minimum radius of 

curvature), 

Z{x'{s),x'{t)) < f \x"(u)\ du < ^— ^ . 

J s r 

Since < s < t < rrr, and the cosine function is decreasing on [0, tt], we have 

(t-s) 

cos(Z(x'(s), x'(t))) > cos . 



Thus 

rt (t — s ) 

(x(t)-x(0))-(x'(t))> / cos^ ds = r shift /r) . 

Jo r 

For < t < m, sin(t/r) > 0. 

Lemma 8 Let K be a C 2 smooth curve in M. 3 , with minimum radius of curva- 
ture r. Let x : [0,£(K)] — > M 3 be a unit speed parameterization of K. Suppose 
0<a<b<c<d< nr, so x(a) , x(b) , x(c) , x(d) are four points in order along 
K , contained in an arc of total length < nr. 
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Then the minimum spatial distance between line segments (x(a)x(b)) and (x(c)x(d)) 
is realized at the closest endpoints. Taking into account Lemma 7, this says, 



MD{ (x(a)x(b)) , (x(c)x(d)) ) = \x(c) - x(b)\ . 



PROOF. Without loss of generality, rescale the curve to have r — 1. Then 
the four points lie in an arc of total length < ir. 

Let A denote the segment (x(a)x(b)) and C the segment (x(c)x(d)). We shall 
show that for each point x G A, the point x(c) is the closest point of C to x; 
so x(c) is the closest point of C to A. By a symmetric argument, the point 
x(b) is the closest point of A to C. 

Fix a point y G C,y ^ x(c),x(d). For any x G A, construct the directed line 
segment from x to y. We claim that the vectors satisfy 

(y — x) • (x(d) — x(c)) > . 

If this dot product is positive, then moving y along C closer to x(d) will in- 
crease the distance to x, and moving y closer to x(c) will decrease the distance 
to x. Thus x(c) must be the closest point of C to x. 

We now show that the above dot product is positive for each x,y. It is con- 
venient to think for a moment of fixing y and varying x. Let P y be the plane 
through y perpendicular to C. Rotate the entire ensemble so that the vector 
x(d) — x(c) points "up". Then the dot product inequality is equivalent to the 
assertion that the entire line segment A lies below P y . It suffices to show that 
each vertex x(a),x(b) lies below P y . 

But in fact, if x(a) and x(b) lie below P x ( c ), then they lie below P y . We have 
now reduced the lemma to the following claim, an inequality that involves 
only the given points on K. The inequality is stated for parameter value a, 
and is identical for b. If < a < c < d < ir, then 

(x(c) - x(a)) ■ (x(d) - x(c)) > . 

The rest of the proof is similar to the proof of Lemma 7 with some trigonometry 
at the end. We first express the difference vectors as integrals of derivatives, 

(x(c) - x(a)) ■ (x(d) - x(c)) = f C [ d x'(s) ■ x'(t) dt ds . 

J a J c 



9 



Since the cosine function is decreasing on [0, 7r], we have 



cos(Z(x'(s), x'(t))) > cos(t — s) . 



As in the proof of Lemma 7, x f (s) ■ x'(t) > cos(t — s), so 



J a 



•c 



■d 



(x(c) — x(a)) ■ (x(d) — x(c)) > 



cos(t — s) dt ds . 



The integral evaluates to 



cos(g? — c) — cos(d — a) — 1 + cos(c — a) , 



which is positive. 

The previous two lemmas tell us that for arcs that are near each other in 
arc-length along a curve, the minimum spatial distance between the arcs is 
the same as the minimum distance between their inscribed chords. For more 
general pairs of arcs, the minimum distances usually will not be equal, but 
they still are related. 

Lemma 9 Suppose a, (3 are smooth arcs in R 3 ; each of length 5, and each 
having radius of curvature everywhere > r > 5. Let e be the chord joining the 
endpoints of a and f the corresponding chord for f3. 

(a) The maximum distance between a and e (likewise between (5 and f) is 

— V48 r 




(b) IfMD(a, f3) is the minimum spatial distance between a and 13, and MD(e, /) 
is the minimum distance between the chords, then 



PROOF. For part (a), imagine the chord e as a rod with "string" of length 
5 attached at either end, and ask, "What configuration allows the string to 
reach as far as possible from the rod?" The answer is when the string is pulled 
out to form two equal sides of an isosceles triangle, with the rod as the base. 
The maximum distance that any point of a can be from e is the altitude h of 

this isosceles triangle, so h 2 = (0 — (y) . Since 5 < r, in particular 5 < nr, 
we can apply Schur's theorem: By Lemma 6 and Lemma 3(b), 




MD(e,/)-MD(a,/3)| < 



< r . 

6r 6 



4 1 1 - 4 48 r 2 



2 ' 
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so 




Part (b) follows from part (a), the triangle inequality, and the fact that 5 < r. 
3. 3 Lemmas about the thickness of a curve 

The first lemma is a characterization of the thickness radius r(K) in terms of 
curvature and the critical self- distance. 

Fix a point x G K and consider points y that start at x and gradually move 
along K. A point y is a critical point for the function \y — x | when y = xq or 
when (xy) _L y'. We define the critical self- distance of K (an idea attributed 
by J. O'Hara to N. Kuiper) to be 



Lemma 10 The thickness of a smooth knot is bounded by the minimum radius 
of curvature and half the critical self-distance. In fact, 



The next lemma is a consequence of Lemmas 6 and 10, and is proven in (1 



Lemma 11 Suppose K is a smooth knot of thickness radius r(K) = r. For 
any x,y e K with arc(x, y) > irr, we must have \y — x\ > 2r. 

Lemma 12 Let K be a C 2 smooth closed curve in IR 3 , with minimum radius 
of curvature r. Let C be a circle whose total arc-length is the same as K , and 
R be the radius of C. Then r < R and (from Lemma 10) the thickness radius 
r{K) < R. 

PROOF. Since r is the minimum radius of curvature of K, the maximum 
curvature of K is -, so the total curvature of K is at most On the other 
hand, by Fenchel's theorem (fljf). the total curvature of K is at least 2tc. Thus 
2vrr < £(K) = 2ttR. 



sd(K) 



min { \y — x 



x 7^ y G K and (xy) _L y } . 




PROOF. See 0). 




11 



K 



P 




C Q 



Fig. 4. The objects of study: smooth knot K with arc «j and vertex Vi, inscribed 
polygon P with vertex Vi, circle C with arc f3i and vertex bi corresponding to on 
and Vi respectively, and inscribed regular ra-gon Q with vertex hi. 

Lemma 13 For any C 2 smooth closed curve K, E L (K) > 2n. 

PROOF. By Lemma 10, the curvature of K is everywhere < l/r(K). Thus 
the total curvature of K is < £(K)/r(K) = E L (K). But the total curvature 
of a closed curve is > 2n. 



4 Notation and outline of proof of Theorem 1 

We have four objects of interest: the knot K, the circle C, the inscribed n-gon 
P, and the regular n-gon Q. In the following list, refer to Figure 4 

• K is a C 2 smooth simple closed curve in M 3 . 

• C is a circle with total arc-length £(C) = £(K). 

• r(K) is the thickness radius of K. 

• K is subdivided into n arcs of equal length 6 = and we are assuming 
5 < r(K) (so n > E L (K)). 

• Vi, . . . , v n are the subdivision points along K. 

• ojj is the arc of K with endpoints Vi and t> i+1 . We number the vertices modulo 
n, so a n is the arc from v n to v±. 

• R is the radius of C, so R = ^P-. 

• t — > x(t) is a unit speed parameterization of K from C. 

• &!,...,&„ are evenly spaced points along C such that = i> «. 
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• Pi is the arc of C corresponding to ctj. 

• P is the polygon formed by connecting the points i>j in order. 

• ti is the edge of P from vi to v i+ ±, with length denoted \ei\. 

• Q is the regular polygon inscribed in C, with vertices bi, . . . ,b n . 

• is the edge of Q with vertices &«, bi+i, with length 

Just to have all the important parameters specified in one place, we also define 
two integers, m and p, whose role will be evident later in this section. 

• m = ^ 7rr C K ") j _ p or a ver tex the vertices i>j,i>j+i, . . . , Vi +m are a maximal 
list that lie in an arc of K of length < nr(K). 

• p = \m*\. For a list of m vertices as specified in the previous item, we 
will need to distinguish an initial bunch from the rest. It turns out that 
the number we need to separate off should be some fractional power of m 
strictly greater than 1/2, and we take 3/4 for simplicity. 

We shall analyze the energies in terms of individual pairs of arcs and/or edges. 

The energies are 

E (K) 



xeK Jy&K \x — y\ 2 \s — t\ 

n n 
i=l j=l 



where 



and 



E (a i ,a j ) = I 



y&aj \X — y\ 



12 ' 



E md (P) = u' md {P) - u'UQ) 

n n 

= E E U md (e u e 3 ) - Umd(fi, fj) (j ^ i - + 1) . (2) 

i=l j=l 

Sometimes we need to treat E as the difference between two integrals, so we 
also define 

Efa, a 3 ) = / 

J x£on Jy^cij \X — y\ 

and likewise for E(/3i, /3j) for arcs on C. As one might expect, our overall plan 
is to show that the various terms in the sum (1) are close to the corresponding 
terms in (2). However, some terms in (1) have no corresponding term in (2); 
and even when they do, there are different cases requiring different analyses. 
We shall, in fact, consider four kinds of pairs bound each contribution 

to the error, and add them to get a full error bound. 

Here is a "schematic diagram" of our situation: We want to show that some- 
thing of the form / (W — X) is close to something of the form (Y — Z). For the 
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edge pairs where Eq has a contribution and E md is not defined, we show the Eq 
contribution is small. For other edge pairs, we sometimes show that J(W — X) 
and (Y — Z) each is small, and sometimes show that \Y — J W\ and \Z — J X\ 
both are small. The analysis has to involve this kind of complication because 
the unregularized polygon energy U' md (P) is not a good approximation of the 
divergent integral JJ KxK \ y l x \i ? that is \Y — J W\ does not get negligibly small 
for arc pairs (and their corresponding segment pairs) that are extremely close 
together along K. Here is a simple example to illustrate the difficulty: Con- 
sider two segments A = [0, e] and B = [2e, 3e] C K. Then U md (A,B) = 1. 
On the other hand, J xeA J yeB , }, 2 dy dx = In | . For segments close together 



along the curves, we need to understand the regularizing terms rather than 
show the two energies are close to each other. 

Following are the four types of pairs (of indices edges or arcs) that de- 

termine our four "zones" for separate analysis. The definitions are symmetric, 
so and are of the same type. 

(1) Adjacent Pairs: j — i — + 1 

For these arc pairs, we bound J2i,j E (ai,aj). Since U md is only defined 
for non-adjacent edges, there are no corresponding edge pairs for these 
arc pairs. 

(2) Near Pairs: non-diagonal pairs for which the arcs a, and otj are 
contained in an arc of K of length < nr(K). 

Within the Near Zone, we make an additional distinction between "Very 
Near" and "Moderately Near": For each vertex Vi, let A be either of 
the arcs of K starting at Vi and having length £(A) = nr(K). The ver- 
tices contained in the arc A are a sequence i>j, i>j+i, . . . , Vi +m (for the 
other arc, we count in the other direction). The arcs contained in A are 
ctj, . . . , aj +m _i. The vertex Vi +m may or may not be an endpoint of A. 
We distinguish between the first m 3//4 vertices and the rest. 

A. For j — i + 2, . . . , i + p, we call (i, j) a very near pair. 

For such we bound (U md {ei, eA - U md {f u fj)) and E {a u aA. 

B. For j — % +p + 1, . . . , % + m — 1, we call (i, j) a moderately near pair. 

For such we shall bound J2ij (E(a h <x,) - U md (e h e^)) and J2i,j (^(A, Pj) ~ U md (f i: }))) . 

(3) Far Pairs: The pairs (i, j) that are neither adjacent nor near are called far. 

For such pairs, we shall also bound E(a h aij) - U md (ei, ej) and £ij E(f3i, fy) - U md (fi, fj), 
but we need an argument different from the moderately near pairs. 

See Figure 5 for an example of the zone pairings where m = 17. We use the 
same terminology for corresponding pairs of arcs in C; that is, if are far 
[resp. adjacent, very near, moderately near] on K, then we call them far [resp. 
adjacent, very near, moderately near] on C. 

In the next section, we establish the explicit error bounds in each of the 
different zones. In Section 6, we collect all of the errors to determine the total 
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Far Zone starts here 



Fig. 5. The four types of zones on which we do our analysis. Note that this is just 
a schematic to show the arrangement of the zones with respect to a fixed arc. 

error bound. 



5 Proofs for the different zones 



5. 1 Bounds for E in Adjacent and Very Near Zones 



We establish the error bound for the combined contributions of the Adjacent 
and Very Near Zones to the Mobius Energy. 

Proposition 14 In the Adjacent and Very Near Zone, 

i,3 



< 1.06 



n 



1/4 



PROOF. If x, y are contained in diagonal or very near arcs, then arc(x, y) < 
(p + 1)5. Thus it suffices to bound 



x&K Jy&K, arc(x,y)<(p+l)5 \X — y\ 



s-t 



The calculation is independent of the choice of x, so we analyze 



\1{K) f 

Jy 



x+(p+l)<5 



\x-y\ 



s-t 



dy 



(3) 



where the limits of integration are meant to indicate that we are integrating 
along an arc of K of length {p + 1)5 starting from x. 
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We are going to find upper and lower bounds for the integrand, observe that 
the upper bound is positive and the lower bound is negative, and conclude 
that the magnitude of the integrand is bounded by the difference between 
the upper and lower bounds. To simplify subsequent expressions, let r denote 
r(K) and a denote arc(x, y). 

Since s and t lie on a circle of radius R, 

1 1 



\s-t\ 2 R 2 (2-2cos(a/R)) 



First we get the upper bound. Since 8 < r, in particular m > 2, we have 
p < m and {p + 1)5 < nr. Thus we can apply Lemma 6 to conclude 

\x-y\ 2 > r 2 (2 - 2cos(a/r)) . 



So we have 

1 1 



< 



\x-y\ 2 \s-t\ 2 ~ r 2 (2-2cos(a/r)) R 2 {2 - 2 cos(a/R)) ' 

By Lemma 12, r < R. By Lemma 6 applied to circles of different radii, or the 
argument in Lemma 5, this upper bound is nonnegative. 

Now we get the lower bound. Since arc-length on any curve must be at least 
as large as chord length, 



Thus 



I |2 ^ 2 

\x — y\ < a . 



1 1 

> — - 



\x — y\ 2 \s — t\ 2 a 2 \s — t\ 2 

1 1 



a 2 i? 2 (2-2cos(a/i?)) ' 
which is negative since chord length < arc-length on a circle. 

Taking the difference between the nonnegative upper bound and the negative 
lower bound, we have 



1 



x — y\ 2 \s — t\ 



1 1 

— r 2 (2 — 2 cos(a/r)) a 2 . 



So 

r(p+i)S 1 1 
(3)<2£(K) / — — — -——-—da, 

K J ~ K J Jo r 2 2 - 2 cos a/r a 2 ' 
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where now we are just integrating a function of a real variable. Applying 
Lemma 4(b), we have 

(3)<2£(K)(p+l)5 1 ' ! ' 1 



4 7T 2 / r 2 

Since m > 3, p > 2, so (p + 1) < 1.5p. Combining the constants, we have 
(3) < 0.45 < 0.45 < 1.06 

V / r 2 - r 2 - n l/4 

as desired. 

5.2 Bound for E rn( i in the Very Near Zone 

Proposition 15 In the Very Near Zone, 

\EUP)\ = PUP) ~ U' md (Q)\ < 2.76 . 



PROOF. 



n t+p 



E md ( yeT J near ) 



\fi\ \fj 



2 ^ 1 ^ +2 MD(e i !e j y~ MD(f\J)) 2 



We shall bound the inner sums uniformly in i, that is bound 



p-i 



E 



Fi| \e i+ k+i\ 



\fi\\f; 



i+k+l\ 



\ MD(e i} e i+k+1 ) 2 MD{f i J i+k+1 f 



(4) 



for arbitrary i, and then multiply that bound by 2n. Here k = j — i — 1 is 
the number of edges separating the two edges. As in Proposition 14, we find a 
positive upper bound for each difference term, and a negative lower bound; so 
the difference between the upper and lower bounds is a bound for the absolute 
value. 



On the circle C of radius R, the edge lengths are \fi 
and MD(fiJ i+k+1 ) = JR 2 (2 - 2 cos(kS/R)). So 



\fi 



R 2 {2-2cos{S/R)) , 



(4) 



p-i 

E 



\ e i\ | e i+fc+l| 



R 2 {2-2cos{S/R)) 



\ MD(e i: e i+k+1 ) 2 R 2 {2 - 2 cos(k5/R)) 
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To simplify subsequent expressions, let r denote r(K). If we compare K locally 
with a circle of radius r, Lemma 8 and Lemma 6 say MD{e^ ei + k+i ) 2 > r 2 (2- 
2 cos(kS/r)). The longest an edge can be is the arc-length, so (|ej| \ej\) < 5 2 . 
Thus, an upper bound for each summand is 

5 2 R 2 (2-2cos(5/R)) 

summand < 



r 2 (2 - 2 cos(M/r)) R 2 {2-2 cos(kS/R)) ' 

We claim this upper bound is positive. First, 5 2 > R 2 (2 — 2cos(5/R)) since 
arc-length (now on the big circle C) is always > chord length. Furthermore, 
r 2 (2 - 2cos(M/r)) < R 2 {2 -2cos(kS/R)) by Lemma 5. 

We next obtain a lower bound. By Lemma 10, r < minimum radius of cur- 
vature of K. So we can apply Lemma 6 and Lemma 8 to any points that lie 
in an arc of K of length < nr. By Lemma 6, we have (|ej| |e i+ fc + i|) > r 2 (2 — 
2cos(5/r)). For the denominator, Lemma 8 gives us that MD(ei, e i+ j, + i) = 
\vi + k+i — Vi + i\, the distance between points of K whose arc-distance is kS. 
Since chord length < arc-length, we thus have MD(ei, ei + t+i) 2 < (k5) 2 . So a 
lower bound for the summand is 

, r 2 (2 - 2 coa(8/r)) R 2 (2-2 cos(5/R)) 
summand > — 



k 2 5 2 R 2 {2-2cos{k5/R))' 

Comparing numerators and denominators as we did for the upper bound, we 
see that this lower bound is always negative. 

Thus, we can bound the absolute value of the summand by the difference 
between the upper and lower bounds: 

5 2 r 2 (2-2cos(5/r)) 

summand < 



r 2 (2-2cos(k5/r)) k 2 5 2 
1 / k 2 5 2 _ r 2 (2-2cos(5/r))~ 

P l v r 2 (2-2cos(A;(5/r)) I 2 ( 



(5) 



We now appeal to our lemmas on cosines and chords. To clarify how lemmas 
will be used, introduce angles 6 = S/r and <fi = k5/r. Thus, the bound (5) can 
be written 

F l v 2-2cos0 I 2 J ' 

By Lemma 3(d), 2 _f COS(j> < 1 + 10 2 . By Lemma 3(c), 2 ~ 2 ™ se > 1 - ±6 2 . Thus, 

Isummandl < — ( -<h 2 + — 9 2 
1 1 ~ k 2 \2 V 12 
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We return to the original double sum and see that 



p-i 



Gi\ |ej+fc+i| \fi\ \fi + k+i 




2n Y, 



MD{ 



Y MD(fi, fi+k+i) 



fc=i 




) 



7 P 




n l/4 



Remark. For the Very Near Zone, we could use p < any fractional power m q . 
It is in the Moderately Near Zone that we need p > 1/2. 

5.3 Bound for \E (K) — E md (P)\ in the Moderately Near Zone 

In this section, we determine the error bounds in the Moderately Near Zone 
for \E(K)- U' md (K)\ and \E(C) - U' md (C)\. Recall that the Moderately Near 
Zone consists of pairs where a, , otj [resp. $ , f3j] are contained in an arc 
of K [resp. C] of length irr(K) but are separated by at least p other arcs; that 
is k = j — i — 1 runs from p to (m — 2). The keys to the analysis in this zone 
are: 

• The minimum distance between a given pair of arcs, or a given pair of 
chords, is realized at the closest endpoints along the curve. 

• That vertex-to-vertex distance is bounded away from zero by Schur's theo- 
rem. 

Proposition 16 In the Moderately Near Zone, 



PROOF. As before, we use r to abbreviate r(K). We first analyze the error 
on K, 



total errorl < 3.00 



n 7/4 



+ 542.84 



E L (Kf/ 2 

n l/2 
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Note: The expressions seem more clear if we use both k and j, where j = 
i + k + 1. 



As in the previous case, the analysis is independent of i, so we work with a 
general i and multiply that bound by n. To bound the above sum of differences, 
we introduce a third term (larger than each of the two we are studying) and 
use the triangle inequality. 

Claim 1. 



m-2 
k=p 



5 2 



MD( ei , ej ) 2 MD(e t , ej y 



< 1.50 



E L {K) 



11/4 



n 



7/4 



Claim 2. 

m-2 
k=p 



MD(e h ej) 2 Jxeat Jy£ aj \x — y 



dy dx 



< 271.42 



E L {Kfl 2 



n 



1/2 



Proof of Claim 1. Since chord length < arc- length, |ej||ej| < 5 2 . So the sum- 
mand without absolute value is non-negative, and any upper bound will bound 
the absolute value. 

Since we are still within the Near Zone, Lemma 6 and Lemma 3(b) give 

1 5 4 



\eil\ej\ >r i (2-2cos(<5/r)) > 5 2 - 



12 r 2 



Now consider the denominator. By Lemmas 8, 6, and 3(b) 

MD(e h ej) 2 = \ Vj - v i+1 \ 2 > r 2 (2 - 2 cos{k5/r)) 



> k 2 5 2 - 



1 k 4 5 4 



Thus 



m-2 / p 

271 £ \MD(e t ,e,) 2 



< - n5 



m-2 
2 \ ^ 



MD{e i ,e j )' 1 J ~ 6 



£ 2 Li-2A2 



k=p 



12 



k 2 5 2 ) 



We next bound this denominator away from 0. In the Near Zone, kS < irr, so 
r 2 — jok 2 S 2 > r 2 (l — ^7T 2 ), which gives 



1 1 1 
< 0.94 



6r 2 (l- jW 2 



2 ' 
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We thus have 



™~ 2 ( S 2 \ ei \ \ ej \ \ 

2n t v \MD(ei, ej) 2 MD^etf) ~ ° MU r 2 ^ F 

5 2 00 1 

k=p 

5 2 1 

< 0.94n — . 

r A p — 1 



We want to bound ^-j- in terms of 5 and r. Recall that p = [m 3 ^ 4 \ and 
m = L^J- Since 5 < r, we have m > 3 and p > 2. So ^-j- < and 
— < Thus, 

m — 3 m+1 ' 

JL < _L. < 3 < 3 /iV" 1 < 3J3 (i_f < L59 f^f . 
p-l-p+1 m 3 / 4 ~ V3/ (m + 1) 3 / 4 \nr J \r J 



Then 

m— 2 



2 - E ( 777v^2 - TTtE^) < (0.94)(1.59)4 



1.50 ^S^. 



< — n 7/4 



This completes the proof of Claim 1. 

Proof of Claim 2. 
We need to bound 



m-2 

fc=p 



MD(e i7 ej) 2 Ja x Jaj \x — y 



12 



(6) 



By Lemma 8, MD(ei, Cj) = MD(oii, aij) = \vj — Since the arcs have 

length 5, we know that the summands without absolute value are nonnegative; 
so, as in Claim 1, we bound the absolute value by finding an upper bound. 
We are dealing with something that looks like a Riemann Sum upper estimate 
of a finite integral. But as n increases, we are changing the domain, not just 
subdividing the same set and we want to control the size of the error, not 
just say it goes to zero as n — > 00. This is where we use the choice of p as a 
fractional power m q where q is strictly larger than 1/2. 
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For brevity, let md denote MD(oti, otj) = \vj — where arc(t>j+i, Vj) = k5. 
Since 

\x — y\ < md + 25 , 

we have 

1 , 5 2 5 2 ..3 1 



< - ~. : =777 < 45 



md 2 Ja t J aj \x — y\ 2 md 2 (md + 25) 2 md 3 

Thus 



(6) < 8n5 3 J2 



m—2 -y 



k=p md3 ' 

As before, by Lemmas 8, 6, and 3(b), since axc(v i+ i,Vj) = k5, 
md 2 = \ Vj -v i+1 \ 2 > r 2 (2-2cos(k5/r)) 
> k 2 5 2 1 k 6 



k 2 S 2 1 



12 r 2 
1 k 2 8 2 " 



12 r 2 



Thus, md > 0.42 k6, so 



> k 2 5 2 ( 1 Ti 2 ) since k5 < nr. 

12 



' < 13.42 "' 



md 3 k 3 5 3 
With the above observation, the £ 3 's cancel and we have 

m— 2 j 

(6) < 107.36 n£- 

k= P K 
00 j 

< 107.36n£- 

< 107.36 n 



(p-1) 2 • 

3/4 



1 f H\ 1 

We showed in the proof of the prior claim that ^-j- < 1.59 (^-J . Thus, 

1 EriK) 3 ' 2 

(6) < 107.36 n- - < 271.42 L \ > — 

y ' - p _ 1)2 - n i/ 2 



Note that in the above analysis the exponent 3/4 needs to be strictly greater 
than 1/2, so that when we double it, the power of n in the denominator will 
more than cancel the leading factor n. 



22 



We now need to bound the contribution from C, that is \E(C) — U' md (C)\. The 
radius of C, R, is the thickness radius r(C). Also, we know from Lemma 5 
that R>r. So if arcs a iy aj of K are near, then the corresponding arcs f3j 
lie within an arc of C of length < irR. Thus, the various steps in our analysis 
of K can be carried out on C. We could obtain sharper bounds for C, but we 
will settle for the same bound since they dominate anyway. 

For Claim 1, we have 

\f t \ |/,-| = R 2 (2 -2cos(5/R)) >5 2 - ±5*R* > 5 2 - ~ , 

and 

MD(fi, f,) 2 > k 2 5 2 - — ^ > k 2 5 2 - — , 
\J»jj) - 12 R 2 ~ 12 r 2 

exactly as for K. Now continue the proof of Claim 1 verbatim. 



For Claim 2, 



mdifijj) 2 = R 2 {2-2cos{k5/r)) 

1 k 2 5 2 \ 




and the rest follows verbatim. 

Thus, our final bound for the total error in this zone is just double the values 
obtained in Claims 1 and 2. 



5.4 Bounds for \E (K) — E md (P)\ in the Far Zone 



As before, we use r to abbreviate r(K). In the Near Zones, we just needed 
a value for r < minimum radius of curvature of K. But in the Far Zone, we 
need both aspects of the thickness radius. 

The argument here is somewhat similar to the Moderately Near Zone, but we 
control the denominators in a different way. In each situation, we need to know 
that spatial distances between points are bounded away from zero in some 
way depending on their arc-length distances along K. For the Far Zone, we 
use the fact that thickness controls critical self-distance, in particular Lemma 
11, together with local analysis (Lemma 9), to relate chord-chord distances to 
arc-arc distances. Also, we continue to use the hypothesis 5 < r. 
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Remark on notation heuristics. In the following paragraphs and Lemma 17, 

think of (cx,f3) as (ccj,aij) and (e, /) as (e^e.,). 

Lemma 17 Suppose (a, (3) is a pair of far arcs (on K or on C), with (e, /) 
the inscribed chords joining their endpoints. Then 

md{a,(3) > 1.08 r , 

and 

md(e,f) > 0.79 r . 



PROOF. We analyze K, and note that the same bound will work for C since 
r < R. We establish the lower bound for arcs, then use that to bound the 
distance for chords. If the minimum distance between a pair of arcs is realized 
at points that are interior to one or both arcs, then we are dealing with singly- 
or doubly-critical pairs of points, so, by Lemma 10, md(a,[3) > 2r. Thus we 
just need to bound the end-point distances. Let ceo and a\ be the endpoints 
of the arc a and [3 and f3\ the endpoints of the arc (3. Choose the labels so 
that cci and (3q are the points which are closest with respect to arc-length. In 
the worst case, the arc-length from a to (5\ is > nr, but the arc-lengths of 
the arcs a f3 , ceiPo, and a.\fi\ are less than nr. In such a case, we have the 
following situation: 

• \a - p \ 

it r > &rc(a , (3 ) > nr — 5 =^ \a — (3 \ 2 > r 2 (2 — 2 cos(7r — 1)) by Lemma 
6, and the fact that 5 < r. So \ao — 0o\ > 1.75r. 

• K - Pi\ 

arc(ao,/3i) > nr \a — j3\\ > 2r by Lemma 11. 

• |«i - Pi\ 

same bound as \ao — fio\ ■ 

• \ai - p \ 

arc(a 1)/ 3 ) > nr — 25 arc(ai, / 9o) > (n — 2)r, since 5 < r. Thus, by 

Lemma 6, 

K~A)| 2 >r 2 (2-2cos(7r-2)) =>■ \a ± - f3 \ > l.08r. 

In other scenarios, the arc pair (a, P) yields three of the above four cases, but 
we lose the smallest. For "most" arc pairs (a, (3), we have all point-to-point 
distances at least 2r. 

We now obtain the lower bound on chord-to-chord distances using Lemma 9: 
md(e, /) > md(a, (5)-^-r> 1.08 r - ^^"j r > 0.79 r . 
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Proposition 18 The total error in the Far Zone is bounded by 



0.56^ + 1.60^ + 7.76^0- 



rr 



n 



PROOF. We first analyze the error on K, 



n n 

2E E 

i=l j=i+m 



md(ei,ej) 2 JxeaiJyeaj \x — y 



dydx 



md(ei,ej) 2 



NN to ,/ 2 ct , that to * 2 , 2 , and 



We do this in three steps: Compare 
that to // 7-3-12. After we do each step for K, we double that to include the 

\x y\ ^ 

contribution from C. Note —J 1 — ^ = L c „ f„ c „ —77-^ — yrdydx and similarly 
for 52 



Claim 1: 



2E E 

i=l j=i+m 



5 2 



md(ei, ej) 2 md(ei,ej) 2 



< 0.28 



E L {KY 
n 2 



Claim 2. 



n n 



2£ E . 

i=i j=i+ m J:rea ' J ^ e % 



md(ei,ej) 2 md(ai,aj) 2 



dyd.v < 0.80 / ' / (/ ° : ' . (7) 



Claim 3. 



n n r r 
i=l j=i+ m J;ce ^ J 2/ ea i 



md(ai,aj) 2 \x — y\ 



dydx < 3.88 v ; . (8) 



n 



Proof of Claim 1. 

Since arc-length > chord length, each summand is nonnegative without taking 
the absolute value, so we just need to bound the terms from above. By Lemma 
6 and Lemma 3(b), 5 2 — < \ei\, \ej\. Thus, 

S 2 |ej| \ej\ 1 5 4 

md(ei,ej) 2 md(ei,ej) 2 ~ 12 r 2 md(ei,ej) 2 

But Lemma 17 gives us that md(e i) ej) 2 > (0.79) 2 r 2 , so 



md(ei,ej) 2 md(ei,ej) 2 r 4 
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Multiplying by 2n 2 gives 



2 EE 

i=l j=m 



md(ei,ej) 2 md(ei,ej) 2 



< 0.28 = 0.28^ 

r 4 n 



Proof of Claim 2. 

The sum (7) is bounded by (2n 2 S 2 ) (worst error in integrands). We will use 
Lemma 9(b) to bound that. To make the algebra more evident, let e = 
md(ei, tj) and 7 = md{a^ otj). The term we wish to bound is 



1 1 




7 2 -e 2 


e 2 7 2 




e 2 7 2 



<L38 i7-ei(7 + e) <a4Q (7 + ^ 2 



since e > 0.79r and 7 > 1.08r by Lemma 17, and I7 — e| < by Lemma 

9(b). 

Now e, 7 are minimum distances between sets that include points of K, so 
e, 7 < e(K)/2 and 7 + e < £(K). Thus, 



e 2 7 2 



< 



0.40£(K)5 2 



Multiplying by 2n 8 , we get 



n n 



2 E E , . 

i=l j=j+m J:ce ^ J V Ga J 



md{e i) ej) 2 md{a i: aj)^ 



dydx < 0.80 



0.80 



Proof of Claim 3. 

The sum (8) is bounded by (2n 2 5 2 ) (worst error in integrand). 

Let 7 denote md(ai, otj). So for particular x, y on a; L and a,-, we have \x — y 
7 + 1 for some < t < 25. The largest error is then 

1 1 _ t(2 7 + t) t(2 7 + 1) t(2 7 + 1) 

7 2 (7 + t) 2 ~ 7 2 (7 + t) 2 _ 7 4 < (1.08) 4 r 4 ' 

since 7 > 1.08r by Lemma 17. 
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Now t < 25 and 7 < £(K)/2. Thus, 



t(2 7 + t) 25(£(K) + 28) 



< 



1.08 4 r 4 ~ 1.08 4 r 4 

2 S(nS + 25) 



< 



1.0* 




J.4: 


2 


5 2 


[n + 2) 


1.0* 


I 4 




2 


5 2 


(2tt + 2)n 



1.08 4 r 4 2vr 
5 2 n 



since n > E L (K) > 2n 



< 1.94 



Thus, 



2 E £ . 

i=i j=i+ m J:r e^ •'l/eay 



md(ai,a>j) 2 \x — y\ 



dydx < 3.88 n 2 S 2 



5 2 n 



= 3.i 



E L (KY 



6 Putting it all together 



n 



Here we combine the bounds from the various zones. 
From Propositions 14, 15, 16, and 18, we have 

\E (K) - E md (P)\ < 3.82E L (K) {^^j ' 

+ 3.00 E L (K) (M^) 7/4 

E L (K) \ 

n J 

2 



1/2 



+ 542ME L (K) 
+ 0.56 E L (K) 



2 , E L (K) 



n 



+ 1.60 E L {Kf ( El{K) 



+ 7.76 E L (K) 



n 

,lE L (Ky 



n 



Since E L (K) > 2n > 1, and n > E L (K), we see that certain terms dominate 
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others. So 



E (K) - E md (P)\ <550 



£lW 5/4 

n l/4 



+ 10 



n 



If n > -Ei(K) 11 / 3 , then the total error is less than 560 El ^/\ 
This completes the proof of Theorem 1. 
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